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"O | ABSTRACT 

^ ■ We report the discovery of a very cool brown dwarf, ULAS J003402. 77-005206.7 (ULAS J0034-00), identified in the 
d ' UKIRT Infrared Deep Sky Survey First Data Release. We provide optical, near-infrared, and mid-infrared photometry 
of the source, and two near-infrared spectra. Comparing the spectral energy distribution of ULAS J0034— 00 to that 
of the T8 brown dwarf 2MASS J04151954-0935066 (2MASS J0415-09), the latest-type and coolest well-studied 
brown dwarf to date, with effective temperature T e ff ~ 750 K, we find evidence that ULAS J0034— 00 is significantly 
cooler. First, the measured values of the near-infrared absorption spectral indices imply a later classification, of 
T8.5. Second, the H— [4.49] colour provides an empirical estimate of the effective temperature of 540 < T e g < 660 K 
(±2tr range). Third, the J- and iJ-band peaks are somewhat narrower in ULAS J0034— 00, and detailed comparison 
against spectral models calibrated to 2MASS J0415— 09 yields an estimated temperature lower by 60 < AT c ff < 120 K 
relative to 2MASS J0415— 09 i.e. 630 < T c s < 690 K (±2ct), and lower gravity or higher metallicity according to 
the degenerate combination —0.5 < A(logg — 2[m/H]) < —0.25 (±2a). Combining these estimates, and considering 
systematics, it is likely the temperature lies in the range 600 < T c s < 700 K. Measurement of the parallax will allow 
an additional check of the inferred low temperature. Despite the low inferred T e g we find no evidence for strong 
absorption by NH3 over the wavelength range 1.51-1.56 /im. Evolutionary models imply that the mass and age are in 
the ranges 15-36 Mj up and 0.5-8 Gyr, respectively. The measured proper motion, of (0'.'37±0'.'07)/yr, combined with 
the photometrically estimated distance of 14-22 pc, implies a tangential velocity of ~30 kms -1 , a value consistent 
with expectation for the inferred age. ULAS J0034-00 is significantly bluer than 2MASS J0415-09 in Y-J, so 
future searches should allow for the possibility that cooler T dwarfs are bluer still. 
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1 INTRODUCTION 

In the decade since the first discovery of a T dwarf 
(G1229, iNakaiima et al. I Il995h over one hundred have 
been identified. Most have been discovered in the field, 
in larg e surveys, especially the Sloan Digita l Sky Survey 
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2007a); and the Second Data Release 



m large surveys, especially trie a loan Digital aky aurvey 
(SDSS; ISfaauss et al. Ill999l ; lYork et al. Il2000l) and the Two 
Micron All Sky Surv e y (2MASS; IBurgasser et al. 1 1 19991 ; 
ISkrutskie et al. 1 I2006T ) . IBurgasser et al. I (|2006aT ) provide a 
classification scheme that uses comparison against a set of 
spectral templates over the range TO to T8, and measure- 
ment of a set of spectral indices that quantify the depth 
of H2O and CH4 absorption features in the J, H and K 
bands. The known T dwarfs have temperature s in the range 
700 <: T eff <; 1500 K ijGolimowski et al. Iliooih . To the end 
of 2006 only six dwarfs classified T7.5 or T8 had appeared 
in the literature. ISaumon et al. I ((2006) and ISaumon et al.l 
(2007) made a detailed study of three of these cool T dwarfs, 
combining near- and mid-infrared photometry and spec- 
troscopy to infer effective temperatures of 725 K ^ T e g ^ 
950 K, surface gravities of 4.8 ^ log[g/(l cms -2 )] ^ 5.4, 
and masses of 25 Mj up ^ M ^ 65 A'/jup. 

The temperature range T c h <^ 700 K rem ains unexplored 
at pr esent, but models such as those of iBurrows et al. I 
(2003) can be used to design searches for brown dwarfs be- 
yond T8. Such very cool brown dwarfs are expected to have 
very low luminosities, to be extremely red in z — J, and, 
at least for T c « ^ 400K, to be blue in J — K. Therefore 
they will be faint in both the z and K bands. Shortward 
of 1 jxm, optical searches such as SDSS are clearly unsuited 
for detecting such very cool brown dwarfs. However, a JHK 
survey, such as 2MASS, is also far from ideal, as cool T 
dwarfs will typically be undetected in K. Samples of candi- 
date cool T dwarfs selected as blue in J — H, but detected 
in only these two bands can suff er from a high degr ee of 
contami nation. As described by iHewett et al. I (120061). th e 
Y band (|Hillenbrand et al. 1120021 ; IWarren fc Hewett II2002T) . 
between z and J, is a valuable addition to the J and H bands 
in a s earch for very cool T dwarfs . Fig. 2 of IHewett et al. 1 
(|2006h . and Fig. 1 of lLodieu et al. I (|2007l ), illustrate the fact 
that T dwarfs in the spectral range T3-T8 are much redder 
in Y — J than main sequence stars of similar J — H colour, 
meaning that a broader range of spectral type of T dwarf 
may be identified than in a simple J — H search, and with 
greatly reduced contamination. Whilst there is some uncer- 
tainty about the near-infrared colours of T8+ dwarfs, they 
should lie in a sparsely populated region of the Y — J vs. 
J — H two-colour diagram, and so the extension of the T 
dwarf sequence to later types should be readily identifiable 
using these colours. 

A wide-field near-infrared survey that is substantially 
more sensitive than 2MASS is needed to find such faint red 
objects, and this was one of several considerations which 
led to the develo pment of the UKIRT Infrared Deep Sky 
Survey (UKIDSS; [Lawrence et al.1l2007l ). UKIDSS is a set 
of five surveys using the UKIRT Wide Field Camera (WF- 
CAM; ICasali et "ai]|2007i ) to obtain imaging of ~20 per cent 
of the sk y in a number of band s selected from the ZYJHK 



set - see 



Lawrence et al.l (|2007T ) for details. UKIDSS began 



in 2005, and is estimated to take seven years to complete. 
Three data rele ases have now ta ken place: the Early Data 
Release (EDR; iDve et al. 1 120061 ); the First Data Release 



(DR1: 
(DR2; 

Of the five UKIDSS surveys, The Large Area Survey 
(LAS) is best suited to searching for very cool brown dwarfs. 
The LAS plans to cover ~4000 deg 2 within the SDSS foot- 
print in the YJHK bands and has an average depth (to 
a signal-to-noise ratio (S/N) of five for point-sources) of 
K ~ 18.2. Critically, this is 2.7 magnitudes deeper than 
2MASS, meaning that the 190 deg 2 covered in YJHK by the 
LAS DR1 fields already probes a volume one fifth that of 
2MASS0 

A search for red point-sources in the UKIDSS 
EDR and DR1 LAS data has alread y yielded nine new 
spectroscopically-co nnrmed T dwarfs (jKendall et alj|2007l ; 
lLodieu et al. 1 120071) over the range T4 to T7.5. Here we 
report the discovery of a very cool brown dwarf, ULAS 
J003402. 77-005206. 7 (hereafter ULAS J0034-00), that ap- 
pears to have a later type, lower effective temperature and 
possibly lower mass than any other known T dwarf. The 
discovery process and follow-up photometric, astrometric, 
and spectroscopic observations are described in Section [5] 
In Section [3] we use a 1 — 2.5 /im spectrum to measure the 
object's spectral type. In Section [4] we examine the spec- 
tral energy distribution (SED) of ULAS J0034-00 over the 
range 1 — 8 /im, using photometry in seven bands, to ob- 
tain an empirical estimate of its temperature. In Section [5] 
the near-infrared spectrum is compared to theoretical spec- 
tral models, providing estimates of the source's fundamental 
physical properties, effective temperature, surface gravity, 
and metallicity, and thereby mass, and age, through the use 
of evolutionary models. In Section [6] we report the results of 
higher-resolution spectroscopy in the H band, to search for 
NH3 absorption. Section [7] provides a brief summary of our 
results. 

The UKIDSS Z YJHK photometry is based on Vega 
l|Hewett et al. I I2OO6). and all magnitudes quoted are in the 
Vega system unless explicitly stated as AB magnitudes (de- 
noted by i(AB), z(AB), etc.). Conversions between AB and 
Vega for the SDSS uqriz and UKIDSS ZYJHK bands are 
given by IHewett et al. I (|2006l ). 



2 OBSERVATIONS AND DATA REDUCTION 

ULAS J0034— 00 was identified as a candidate T dwarf after 
a largely automated search of the UKIDSS DR1 database 
(Section 12. ip . after which follow-up photometry confirmed 
that its colours are consistent with that interpretation (Sec- 
tion l2.2[) . Photometry at mid-infrared wavelengths has been 
obtained (Section l2.3[) and the proper motion has been mea- 
sured using a one-year baseline (Section 12. 4p . Near-infrared 
spectroscopy reveals ULAS J0034— 00 to be a brown dwarf 
of very-late spectral type (Section 12. 5[) . 



2.1 Identification 

ULAS J0034— 00 was discovered as a product of a gen- 
eral search of the UKIDSS DR1 catalogue for objects 

1 The gain in depth in the H band, which sets the volume in a 
search for very cool brown dwarfs, is very similar to the gain in 
K. 
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Figure 1. The UKIDSS DM F-band image of ULAS 
J003402.77-005206.7. The exposure time is 40s. The field of view 
is 3' on a side. 

with extreme optical-to-near-infrared colours, including very 
cool brown dwarfs, as well as quasars of redshift z > 6 
ijHewett et al. 1120061 ). The selection was made by a SQL 
query of the WFCAM Science Archive (WSA). A detailed 
description of the W SA, together w i th ex amples of SQL 
queries, is provided bv lHamblv et al~l ([2007) . The first step 
in the search was to match the UKIDSS DR1 catalogue 
to SDSS DR5 (|Adelman-McCarthv et al. Il2007r ), retaining 
both unambiguous matches (within a match radius of 2'.'0) 
with i — Y ^ 3, as well as UKIDSS sources with no SDSS 
counterpart. The sample was further restricted to point-like 
sources with Y ^ 19.5. The sample was cleaned of cross- 
artefacts by rejecting sources located up to seven inte- 
ger multiples of 51'.'2 in RA or dec from bright (J ^ 13.5) 
2MASS stars. 

The UKIDSS DR1 includes LAS fields covering some 
100 deg 2 within a section of the multiply-scanned SDSS 
Stripe 82 (the ~212.5deg 2 define d by -25° < q(2000) < 
+60°, -1.25° < 5(2000) < +1.25°. [Stoughton et al.ll2002j) . 
Therefore it was possible to obtain improved i- and z-band 
p hotometry of candidate s within these fields (as described 
in lVenemans et al. II2007T ). Most of the candidates were re- 
vealed to have colours consistent with their being M dwarfs 
and were thus rejected. 

The resultant list of red point-sources could include 
both high-redshift quasars and brown dwarfs. We used the 



2 Saturated images of bright stars observed by WFCAM pro- 
duce electronic ghost images in nearby detector channels, at in- 
teger multiples of the detector channel spacing of 128 pixels. The 
phenomenon, known as cr o ss-tal k, is illustrated and described in 
more detail in iDve et al. Because ghosts appear at the 

same location in each of the UKIDSS filters they are often cat- 
alogued as being detected in multiple bands, but, since they are 
not real, they have no SDSS counterparts, thus masquerading as 
potentially interesting red sources. 



synthetic colours of He wett et al. I (|2006h to divide the list 
into a sample of candidate quasars, with Y — J ^ 0.8, and 
a sample of candidate brown dwarfs, with Y — J > 0.8. 
Quasars of redshifts 6 < z < 7 have predicted colours 
Y-J~ 0.5, J-H~ 0.4, while mid to late T dwarfs 
have predicted colou rs Y — J~ 1.1, —0.5 < J—H < 0.3 
jHewett et al. 1 120061 '). ULAS J0034-00 was located in the 
sample of quasar candidates, and stood out as the only can- 
didate which remained undetected in any band in the co- 
added SDSS images. 

The Y-band discovery image of ULAS J0034-00 is 
shown in Fig. [T] The UKIDSS DR 1 point-source photom- 
etry (apermag3, IDve et al. I 120061 ) is Y = 18.90 ± 0.10, 
J = 18.14 ± 0.08, H = 18.50 ± 0.22, K > 17.94 (5a de- 
tection limit, established from the measured sky noise). The 
measured colours Y-J= 0.76 ± 0.13, J-H= -0.36 ± 0.23, 
were ambiguous, the Y — J colour suggesting a quasar, and 
the J — H colour suggesting a brown dwarf, and so a number 
of follow-up observations were made. 



2.2 Optical and near-infrared photometry 

To clarify the nature of the source, higher S/N images in 
the z, J, H, and K bands were obtained. The best photom- 
etry for the source across the wavelength range z to K is 
summarised in Table [1] Some additional photometry, subse- 
quently superseded, is provided in the text, below. 

The source was observed for 1800 s with filter #611 in 
EMMI, at the NTT, on the night beginning 2006 August 
18. This filter/detector combination is referred to here as z„ 
and is similar to the SDSS z band. To calibrate the image 
we first used the measured CCD sensitivity curve and the 
filter transmission curve to establish the colour relation be- 
tween z n and SDSS i and z. Briefly, the AB magnitudes in 
the i, z, and z n bands of dw arf stars in the Gunn Stryker 
atlas (|Gunn fc Stryker Hl983T l were computed by appropri- 
ate integration under the bandpasses, to establish the colour 
relation z„(AB) = z(AB) - 0.05(i(AB) - z(AB)) for dwarf 
stars from O to M. Photometry in i(AB) and z(AB) of SDSS 
DR5 sources in the field was then converted to z„(AB), and 
the brightness of the source was measured by relative pho- 
tometry using a fixed aperture, and is recorded in Table [1] 
We also computed the colours of some late T dwarfs with 
good spectrophotometry, to determine that ULAS J0034— 00 
would be about 0.2 mag. fainter in z(AB). 

The offset to the Vega system is z„(AB) = z n + 
0.563, computed by inte grating the spectrum of Vega of 
iBohlin fc Gillilandl (|2004h under the V and z n bandpasses, 
adopting V — +0.026, and demanding z n — V i.e., zero 
colour for Vega. 

The source wa s observed with the UFTI instrument 
l|Roche et al. 1 120021 ') on UKIRT in J and H for 450 s each, 
on the night beginning 2006 August 7, and in K for 1620 s, 
on the night beginning 2006 September 2. The data were 
calib rated using the UKIRT Faint Standards l|Leggett et al. I 
2006). The filters are from the Mauna Kea Observatories set 
(MKO pTokunaga et al. l [2002). and the phot ometric system 
is vir tually identical to the WFCAM system l|Leggett et al. I 
2006). The data were reduced using the astronomical imag- 
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Table 1. Optical and near-infrared photometry of the source 
ULAS J0034-00. 



Band 


mag. 


Instrument 


Date 




z„(AB) 


21.91 ± 0.05 


EMMI 


2006-08- 


18 


Y 


18.90 ±0.10 


WFCAM 


2005-10- 


04 


J 


18.15 ± 0.03 


LIRJS 


2007-01- 


03 


11 


18.49 ± 0.04 


LIRIS 


2007-01- 


03 


K 


18.48 ± 0.05 


UFTI 


2006-09- 


02 



ing pipeline ORAC-DR ijCavanaeh et al. I l2003h PI The re- 
sult for K is provided in Table [1] The other measures were 
J = 18.14 ± 0.03, and H = 18.51 ± 0.05. This are not pro- 
vided in the table as they have been superseded (see below) . 

The blue colour J — H — —0.37 ± 0.06, measured with 
UFTI, established that the source is a cool brown dwarf 
rather than a quasar. It was possible to estimate an ap- 
proximate spectral type by comparing the measured colours 
z n (AB) — J — 3.77±0.06, J-H= -0.37±0.06, with the ta- 
ble of synthetic colou rs and spectral types of T dwarfs from 
IHewett et al~l (|2006h . We first convert z n (AB) to z(AB) for 
ULAS J0034— 00 by adding 0.2 (as above), and then convert 
to the Vega sys tem by subtra c ting t he AB offset of 0.533, 
from Table 4 in IHewett et al. I l|2006h . yielding the approxi- 
mate Vega colour z — J = 3.44±0.06. Taking t he colours and 
spectr al types of T dwarfs from Table 10 of IHewett et al. 1 
(2006), we found there is a correlation between each of the 
two colours and spectral type. From the measured colours of 
ULAS J0034— 00, the errors, and the scatter in the relations, 
we estimated an approximate spectral type T7± 1. Accord- 
ingly the object was ranked as a high priority candidate for 
spectroscopy. 

The target was reobser ved at higher S/N in J and H 
with the LIRIS instrument (|Acosta-Pulido et al. 1120031 ) on 
the WHT, on the night of 2007 January 3. Integration times 
were 900 s and 1000 s respectively. These data were cali- 
brated to the UKIDSS DR1 images using nearby bright un- 
saturated stars. We added an error of 0.02 mag. in quadra- 
ture to the random photometric errors, to account for the 
accuracy of the UKIDSS calibration. The results are pro- 
vided in Table[T] 

The colour Y - J = 0.75 ± 0.10 of ULAS J0034-00 
is rather blue in comparison with other T dwarfs found in 
UKIDSS. The average colour of the nine ne w spectroscop- 
ically confi rmed T4 to T7 . 5 dw arfs listed in iKendall et al.l 
l|2007h and lLodieu et al. I (120071 ) is 1.13. For completeness 
we note that the zero point of the UKIDSS Y photometry is 
under review. Preliminary analysis indicates that a constant 
0.09 mag. should be added to the UKIDSS DR1 Y values, 
to plac e the photometry on the Vega system (| Warren et al. I 
l2007af ). Nevertheless, recalibration will not, of course, alter 
the difference in colours quoted above. 



2.3 Mid-infrared photometry 

We obtained IRAC four-channel (3.55, 4.49, 5.73 and 
7.87 /im) photometry of the source on UT date 2006 De- 
cember 26. The data were obtained as part of the Spitzer 
Space Telescope DDT program #281. All four channels have 
256x256-pixel detectors with a pixel size of l'.'2xl'.'2, yield- 
ing a 5!2x5?2 field of view. Two adjacent fields are imaged 
in pairs (channels 1 and 3; channels 2 and 4) using dichroic 
beam splitters. The telescope is then nodded to image a tar- 
get in all four channels We used exposure times of 30 s 
and a 5-position medium-sized (52 pixels) dither pattern re- 
peated ten times, for a total observing time of 57.6 minutes. 

The data were reduced using the post-basic-calibration 
data mosaics generated by version 15 of the IRAC pipeline^ 
The mosaics were flat-fielded and flux-calibrated using 
super-fiats and global primary and secondary standards ob- 
served by Spitzer. 

We performed aperture photometry using an aperture 
with a 2-pixel (or 2'.'4) radius, to minimise the contribution 
from a faint source 3'.' 5 to the south-east. To assess the influ- 
ence of the neighbouring source, we repeated the photometry 
using an aperture with a 4-pixel radius. After applying the 
appropriate aperture corrections, the results for the larger 
aperture are total fluxes 6 to per cent brighter, for channels 
1 through 4, respectively. Therefore contamination of the 
2-pixel photometry by the neighbouring source will be at a 
level substantially below these values. To convert the 2-pixel 
fluxes to total fluxes, we applied aperture corrections as de- 
scribed in Chapter 5 of the IRAC Data Handbook 3 of 1.205, 
1.221, 1.363 and 1.571 to channels 1 through 4, respectively. 
The photometry was converted from milli-Janskys to mag- 
nitudes on the Vega system using the zero-magnitude fluxes 
given in the IRAC Data Handbook (280.9, 179.7, 115.0, and 
64.1 Jy for channels 1 to 4, respectively). Note that IRAC 
observations are reported as a flux density at the nominal 
wavelengths given in Table [2] assuming that the target has 
a flux density /„ oc \jv. This assumption is not valid for 
T dwarfs and so the results given in the Table should not 
be used to derive a spectral flux at the nominal wavelength. 
Howev er, if the mid-in f rared spectral energy distribution is 
known , ICushing et al. 1 (|2006h show how the IRAC fluxes can 
be used to photometrically calibrate the spectrum - or com- 
pute the 'color correction', in the terminology of the IRAC 
Data Handbook. 

Photometric random errors were derived from the un- 
certainty images that are provided with the post-basic- 
calibration data. The magnitudes and errors are given in 
Table [2] There ar e, in addition, sys tematic errors from a 
variety of sources dReach et alj|200a ). and the IRAC Data 
Handbook recommends citing 5% uncertainty for the abso- 
lute calibration. This includes a contribution from the cal- 
culation of the 'color correction', which therefore does not 
apply here. Additional sources of uncertainty include: the 
accuracy of the ab solute calibration (estimated to be 2% by 
iReach et al.ll2005l ). systematic errors introduced by pipeline 



3 The ORAC-DR Imaging Data Re- 

duction User Guide is available at 

http://www.starlink.ac.uk/star/docs/sun232.htx/sun232.html 



4 For more information about IRAC, see 

iFazio et alH d2004h and the IRAC Users Manual at 
http: / / ssc . spitzer . caltech.edu / irac / descrip.html 
'° Information about the IRAC pipeline and data products can 
be found at |http: / / ssc. spitzer. caltech.e du/irac/dh / 
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Table 2. Mid-infrared (Vega) photometry of the source 
ULAS J0034— 00. A 3 per cent error should be added in quadra- 
ture to the quoted random errors to account for systematics. 



Band 


mag. 


/tin 




3.55 


16.28 ±0.01 


4.49 


14.49 ±0.01 


5.73 


14.82 ±0.04 


7.87 


13.91 ±0.05 



dependencies (of comparable size, iLeggett et al ] |2007al l. the 
variable pixel scale across the field of view, and the variation 
of the quantum efficiency across a pixel. For observations 
confined to near the centre of the array, the last two effects 
combined appear to be at a level of < 2% in an individual 
exposure, as measure d from the scatter in repeat dithered 
observations (Fig. 3 in lPatten et al.|[2006l ). and therefore ex- 
tremely small after averaging dithered observations. Accord- 
ingly, we adopt the total photometric uncertainty to be the 
sum in quadrature of the values given in Table [2] plus 3 per 
cent. 



2.4 Astrometry 

The original DR1 YJHK images were all taken on the 
night beginning 2005 October 4, and serve as a first epoch 
for the measurement of ULAS J0034— 00's proper motion. 
Since the J image has the highest S/N this is the best 
wavelength for a second-epoch image. The source name 
ULAS J003402. 77-005206.7 is derived, following the pro- 
tocol of the International Astronomical Union, from the 
UKIDSS DR1 coordinates measured on the Y image of 
h 34 m 2.771 s - 0°52'6.78". The coordinates are calibrated 
to 2MASS, and are on the International Celestial Reference 
System. The absolute accu racy is better than 100 mas on 
each axis (|Dve et al. 112006] ). 

A deeper WFCAM J image was obtained on the night 
beginning 2006 December 4. The source was observed for 
400 s, and was located on the same detector as for the 
first epoch, in order to minimise systematic errors due to 
distortion. The measured shifts in RA and dec relative to 
the first epoch, and the computed proper motion are pro- 
vided in Table [3] We detect a significant proper motion of 
0"37±0'.'07/yr. 

The LIRIS H-band image (Section 12. 2[) . compared to 
the first epoch H image, provides an independent measure 
of the proper motion. The result, provided in the second line 
of Table [3] is consistent with the J-band measure, but with 
a larger uncertainty. Rather than combine the two results, 
we retain the J-band measure as the best estimate, because 
the differential systematic errors due to distortion should be 
smaller, since the two J-band images were taken with the 
same instrument. 



2.5 Spectroscopy 

2.5.1 GNIRS low-resolution spectrum 

The Gem ini Near-Infrared Spectrograph (GNIRS; 

lElias et al. I [2006) on Gemini South was used to make 



quick response observations of ULAS J0034— 00, through 
program GS-2006B-Q-36. GNIRS was used in cross- 
dispersed mode with the 321/mm grism, the l'.'O slit and 
the short camera, to obtain 0.9-2.5 /im R~500 spectra. 
Two triggers were made on the target, providing a shorter 
and longer total integration time of 16 and 58.7 minutes, on 
the nights of 2006 August 26 and 2006 Oct 1, respectively. 
The target was nodded 3'.'0 along the slit in an "ABBA" 
pattern using individual exposure times of 240 s and 220 s 
respectively. Calibrations were achieved using lamps in the 
on-telescope calibration unit. Early F stars were observed 
as telluric standards, at an airmass that matched the 
mid-point airmass of the target observation. 

Data reduction was initially undertaken with tasks in 
the Gemini GNIRS IRAF Package. The cross-dispersed 
spectrum comprises five highly curved spectra, spatially sep- 
arated, that correspond to different grating orders over the 
wavelength of interest. The reduction procedure included 
correction for offset bias, order separation, median stack- 
ing at the A and B positions, sky subtraction, S-distortion 
correction (i.e. straightening the curved spectra), and wave- 
length calibration. The data were not flat-fielded. Apart 
from bad pixels, the variation of the array quantum effi- 
ciency is sufficiently small that flat-field correction makes 
no discernible reduction in the noise in the spectra. Fur- 
ther reduction was carried out in IDL. Sky residuals were 
fit, and subtracted, using a surface constructed via a se- 
ries of linear fits across the slit (excluding pixels within the 
spectral apertures). Spectra were then extracted within l'.'5 
apertures and summed for the A and B positions. Noise spec- 
tra were also extracted. The spectra were sigma-clipped and 
calibrated onto a relative flux scale using the spectrum of 
the telluric standard (after appropriate interpolation across 
any hydrogen absorption lines) and a black-body function 
for T e fr= 7500 K. The spectral orders were then trimmed of 
their noisiest portions, and the spectrum was scaled to ab- 
solute flux by the J-band photometry. The final spectrum 
consists of an average of the longer and shorter triggered 
observations, weighted by exposure time. 

The pixel scale oversamples relative to critical sampling 
by a factor three (the slit is 6.7 pixels wide). Because the 
spectrum S/N is relatively low, and bearing in mind that 
the principal diagnostic features are broad, for the purposes 
of display the spectrum has been binned by a factor of 12 
(i.e. two resolution elements per pixel). The complete spec- 
trum is shown as the magenta line in the upper panel of 
Fig. [2] where it is compared to the spectra of two well- 
known very c ool brown dwarfs , the T7.5 dwarf Gl 570D, 
discovered by Burgassc r et al. I (l200Ch. and the T8 dwarf 
2MASS J0415-09, discovered by iBurgasser et al. I ll2002ri . 
The comparis o n spe ctra are from lGeballe et al.l (120011 ) and 
iKnapp et al. I (|2004l ). respectively. The comparison spectra 
are of comparable resolution to the GNIRS spectrum. The 
lower panel plots model spectra, and is discussed in Section 
15.21 Fig. 13 provides a zoom of the absorption troughs either 
side of the J and H peaks. 

We synthesised colo urs from the sp e ctrum , following the 
procedure described by iHewett et al.l l|2006l ). Briefly, this 
involves integration of the spectrum under the passband re- 
sponse curves, with appropriate normalisation by the same 
integration performed on the spectrum of Vega. The result- 
ing colours are as follows, and are in satisfactory agreement 
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Table 3. Measured shifts relative to the base epoch of 2005-10-04, and corresponding proper motions. The proper motions (but not the 
shifts) have been corrected for parallax, assuming a distance of 18 pc. 
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epoch 1 


epoch 2 
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tt 


''/yr 


"/yr 


J 


2005-10-04 


2006-12-04 


-0.18 


-0.43 


-0.12 ±0.05 


-0.35 ±0.05 


H 


2005-10-04 


2007-01-03 


-0.09 


-0.39 


-0.04 ±0.07 


-0.31 ±0.07 




wavelength /am 



Figure 2. Upper: GNIRS spectrum of ULAS J0034— 00, binned by a factor 12, plotted magenta. The orange line plots the la error 
spectrum. The five separate sections correspond to different spectral orders. The blue and green lines plot the spectra of the T7.5 brown 
dwarf Gl 570D and the T8 brown dwarf 2MASS J0415— 09, for comparison. These spectra were scaled to the crown of the J-band 
peak of ULAS J0034— 00, over the wavelength range 1.25-1.28 lira. Lower: Solar— metallicity model spectra representing ULAS J0034— 00 
(magenta), Gl 570D (blue), and 2MASS J0415-09 (green). 



A very cool brown dwarf 7 




1.2 1.4 1.6 

wavelength jj,m 



Figure 3. Zoomed portion of Fig. [2]showing the absorption troughs either size of the J and H peaks. Wavelength regions of the spectral 
classification bands (Table [4} are marked. Upper: GNIRS spectrum of ULAS J0034— 00, binned by a factor of 12, plotted magenta. The 
orange line plots the la error spectrum. The separate sections correspond to different spectral orders. The blue and green lines plot 
the spectra of the T7.5 brown dwarf Gl 570D and the T8 brown dwarf 2MASS J0415— 09, for comparison. These spectra were scaled 
to the crown of the J-band peak of ULAS J0034— 00, over the wavelength range 1.25—1.28 /im. Lower: Solar-metallicity model spectra 
representing ULAS J0034-00 (magenta), Gl 570D (blue), and 2MASS J0415-09 (green). 



with the measured values, provided in parentheses: Y — J = 
0.79±0.03 (0.75±0.10), J-H = -0.46±0.01 (-0.34i0.05), 
H - K = 0.05 ± 0.02 (0.01 ± 0.07). As previously noted the 
measured Y — J colour is unusually blue. The synthetic Y — J 
colour corroborates this finding. The synthetic Y — J colours 
of Gl 5 70D and 2MASS J0 415-09 are 0.95 and 1.06 respec- 
tively (|Hewett et al. 1120061 ). 



2.5.2 ISAAC intermediate-resolution spectrum 

An intermediate-resolution spectrum covering the wave- 
length range 1.50 — 1.58^tm was obtained with the ISAAC 
instrument on the ESO VLT using Director's Discretionary 
Time, under programme 278.C-5014(A). The purpose of 
the spectrum was to search for additional absorption lines 
in comparison with the spectrum of 2MASS J0415— 09, 
that could explain the excess absorption seen in the low- 
resolution spectrum of ULAS J0034-00 (Figs [2] and [3| , and 
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might be due to NH3. We used the MR grating and a 1" slit, 
providing a nominal resolving power of ~ 3000. The aver- 
age seeing was 0. 7, meaning that the actual resolving power 
was nearer 4000. The dispersion was 0.807A, corresponding 
to ~ 5 pix per resolution element. 

The source was observed on the nights beginning 2007 
Jan 3 and 11. On each night, integrations of 600s each were 
obtained at four separate slit posit ions, providing a total o f 
80min on source. As described in Wcathcr lev et all (120051 ) 
higher S/N can be achieved with multiple slit positions, com- 
pared to the traditional method of observing in ABBA se- 
quence and subtracting in pairs. For each slit position a sky 
frame was created by averaging the three other frames from 
the same night, and was subtracted. The data were then 
flat-fielded, and improved sky subtraction achieved by fit- 
ting a function up each line. All eight frames were then reg- 
istered to the nearest pixel spatially and spectrally, scaled 
to a common count level, and averaged, ignoring bad pix- 
els. Wavelength calibration was achieved using the sky lines, 
and was found to be linear to better than a pixel. Therefore 
rebinning was not required in applying the wavelength solu- 
tion. Observations of standard stars were used to correct for 
atmospheric absorption (everywhere less than 10 per cent 
over this wavelength region), and to calibrate onto a rela- 
tive flux scale. Absolute calibration was achieved by scaling 
to the GNIRS spectrum. 

The spectrum is discussed later in Section [6] where the 
flux-calibrated spectrum, together with the error spectrum, 
is plotted in Fig. [8] 



3 SPECTRAL TYPE 

To determine the spectral type we follow the revised clas - 
sification scheme for T dwarfs of iBurgasser et al. I (|2006af ) 
which unifies and refines the two independ e nt pr elimi- 
nary classification sch emes of iGeballe et al.l (|2002) and 
IBurgasser et al. I (120021), both of which are grounded in the 
methodology of iMorgan fc Keenan I (|l973l ). Classification is 
based on the depth of the absorption bands due to H2O 
and CH4 in the near-infrared, between the Y, J, H, and K 
peaks, which strengthen towards later spectral types, and 
are nearly saturated by T8. Classification is achieved by 
comparison of the object spectrum against a set of template 
spectra, or from the measured values of a set of spectral 
indices that quantify the strength of the absorption. 

In Figs. [2] and [3] the spectrum of ULAS J0034-00 
is compared against that of the T8 spectral template, 
2MASS J0415-09. Also plotted is the spectrum of the 
well studied bro w n dwa rf Gl 570D, classified as T7.5 by 
IBurgasser et al.l |2006a). The three spectra have been 
scaled to the crown of the J-band peak. There is visibly 
a progressive narrowing of the J-band peak in moving from 
Gl 570D, to 2MASS J0415-09, to ULAS J0034-00, as well 
as successively stronger absorption in the H2O trough near 
1.15 /xm. If the spectra were instead normalised to the H- 
band peak the same sequence would be evident in the blue 
wing of the peak, near 1.5 /im. This visual comparison rules 
out a spectral type as early as T7.5 for ULAS J0034— 00, 
and suggests a classification of T8 or T8.5. 

The K-band peak is relatively stronger in 
both ULAS J0034-00 and Gl 570D, compared to 



2MASS J0415— 09, but this region is not used in cur- 
rent spectral classification schemes. A weak absorption 
feature in the spectrum of ULAS J0034— 00 near 2.0 )im 
coincides with a strong telluric band, and is probably not 
intrinsic to the source. 

Besides comparison against template spectra, we have 
also measured t h e (prim ary) spectral indices defined by 
IBurgasser et al.l l|2006al ). and these are given in Table [1] 
The first three indices imply a classification of T8.5, the 
fourth gives T8 and the fifth, CH4— K, saturates at T7, and 
therefore is not useful for this source. Taken together the 
indices suggest that T8.5 is the appropriate class. This is 
in agreement with the value from the template comparison, 
and therefore is our adopted value. 

Although the precision of the classification scheme itself 
is only 0.5 of a class, this is the first brown dwarf discov- 
ered with a classification later than T8. For this reason we 
have investigated the uncertainty in the class in a formal 
way, treating class as a continuous variable. Then the value 
of each index and its uncertainty, together with the mean 
and range of the parameter defining the T8 class, provide 
four estimates of the parameter. An unweighted mean of the 
four estimates yields a class T8.5±0.1 (error on the mean). 
For the weighted mean the uncertainties on the indices were 
doubled, to force x 2 = 3 (for three degrees of freedom), 
which led to the result T8.4±0.1. A systematic error in the 
sky subtraction, summed over any band, could make a larger 
contribution to the error budget. Nevertheless this analysis 
supports the result from the template comparison that the 
class cannot be earlier than T8. 



4 EFFECTIVE TEMPERATURE ESTIMATE 
FROM THE SPECTRAL ENERGY 
DISTRIBUTION FROM 1-8 /jM 

The SED over the range 1-8 /im, measured by the broad- 
band photometry, provides furt her evidence of the ex- 
treme nature of ULA S J0034-00. iPatten et ail <|2006h and 
iLeggett et al.l |2007al ) note relations between various colours 
within this wavelength range and spectral type, luminos- 
ity, and effective temperature. ULAS J0034— 00 is more 
extreme than 2MASS J0415— 09 in nearly all the colours 
plotted. In Fig. [4] we plot the photometry in the JHK 
bands plus the four IRAC bands for ULAS J0034-00 
and for 2 MASS J0415- 9. Fo r the latter the data were 
taken from IPatten et alJ l|200eh , and we note that, as for 
ULAS J0034-00, the JHK photometry is on the MKO sys- 
tem (no Y-band photometry for this source exists as far 
as we are aware) . The datapoints in Fig. [4] should only be 
interpreted in a relative sense because, as described in Sec- 
tion l2.3l the IRAC magnitudes are not equivalent to the cor- 
responding flux at the nominal filter wavelength. The SED of 
ULAS J0034— 00 is significantly redder over the entire wave- 
length range plotted, indicative of lower temperature. For 
example, colours for ULAS J0034-00 (2MASS J0415-09) 
are as follows; J - K = -0.33 ± 0.06 (-0.51 ± 0.04), 
K - 3.55 = 2.20 ± 0.06 (1.73 ± 0.04), 3.55 - 7.87 = 2.37 ± 
0.07 (1.99 ±0.06). 

The absorption bands of H2O and CH4 which de- 
fine the T spectral sequence are nearly saturated at T8. 
This means that spectral changes in the near-infrared 
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Table 4. Measured values of the spectral classification indices of iBurgasser et aT, (2006a). for the source ULAS J0034— 00. Each index 
is labelled by both the absorption species and the band in which it appears, and is defined by the ratio of the summed flux in the 
two wavelength ranges given. The error associated with the measured index is only the random contribution obtained by propagating 
the uncertainties implied by the noise spectrum. For each index the range associated with the T8 classification is given, along with the 
classification for ULAS J0034— 00 implied by the measured value. 



Band numerator denominator value T8 type 

wavelength wavelength range 
/an fj,m 





-J 


1.14 - 


1.165 


1.26- 
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may be relatively insensitive to temperature for objects 
cooler than 2MASS J0415-09. Therefore it is conceiv- 
able that ULAS J0034— 00 is substantially cooler than 
2MASS J0415— 09, despite the rather small near-infrared 
spectral differences. In order to estimate T e g we have com- 
piled a sample of L and T dwarfs that possess good temper- 
ature estimates, as well as photometry in the seven bands 
from J to [7.87] plotted in FigJU Temperatures a re take n 
from iGolimowski et ai~l (|2004T ); ILuhman et al. I (l2007t): 
iLeggett et alj ll2007bl) . and p h otometry from [Patten et alj 
l|2006t) : ILeggett et all (1200781 ); ILuhman et al. I (120071 ). We 
then plotted T e a against all possible colour combinations, 
looking for the best correlation. We avoided the K-band be- 
cause the flux in this band is sensiti ve to both gravi ty and 
metallicity at fixed temperature (e.g. lLiu et al.| [2007). 

In Fig. [5] we plot T e g against the if— [4.49] colour for L 
and T dwarfs with T c g< 2000 K. The L dwarfs become red- 
der in H— [4.49] towards cooler T e g, but there is a discontinu- 
ity in the relation at the L/T transition, reflecting the strong 
spectral changes that occur over a small temperature range. 
Nevertheless for later T dwarfs with T c ff< 1250 K there is a 
remarkably tight correlation. The RMS deviation from the 
best-fit linear relation, plotted in the figure, is only 28 K. 
The sample plotted includes the peculiar T6 dwarf 2MASS 
J0937±2 9. The K-band flux of this object is strongly sup- 
pressed (|Burgasser et al.~ll2002l). due either to hig h gravity 
or low metallicity or both ( Burgasser et alj|2006bf h The ob- 
ject is the most distant outlier in Fig. [5j located at T e s 
= 900 K, if- [4.49] = 3.00. Nevertheless it lies only 52 K 
above the line, i.e., at 2a. This hints that for mid and late T 
dwarfs the if— [4.49] colour is relatively insensitive to varia- 
tions in gravity and metallicity. 

Other than ULAS J0034-00, the reddest object plot- 
ted in Fig. [5] is 2MASS J0415-09, which has if- [4.49] = 
3.41 ± 0.04 and an effective temperature of T e g — 750 K. 
ULAS J0034-00 is considerably redder, with #-[4.49] = 
4.00 ± 0.05, implying that it is correspondingly cooler. Ex- 
trapolating the linear fit provides a temperature estimate 
T cff = 600 ±30 K for ULAS J0034-00. The lower panel in 
Fig. [5] is a zoom of the region of interest. Also shown are 
model tracks in this parameter space, computed from the 
synthetic spectra described in Section [5] for a range of grav- 
ities (parameterised by logg), and metallicities [m/H]. Over 
the temperature range 700 < Tea < 900 K the model tracks 
fit the data well, while the small scatter between the tracks 
implies that the small scatter seen in the data is not sim- 



ply fortuitous, but because the H— [4.49] colour is indeed 
relatively insensitive to variations in metallicity and grav- 
ity. The model tracks remain linear to cooler temperatures 
T c ft < 700 K, providing some justification for extrapolat- 
ing the linear fit to the data. In fact the tracks are some- 
what steeper than the linear fit to the data, which implies 
a temperature even lower than 600 K. We defer a more de- 
tailed discussion of the models in the mid-infrared region 
until completion of scheduled mid-infrared spectroscopy. 

Assuming the linear extrapolation of the colour- 
temperature relationship to be valid, ULAS J0034— 00 could 
have a comparable temperature to 2MASS J0415— 09 only if 
it were a ~ 5a outlier. In this respect it would be more ex- 
treme than 2MASS J0937±29 (mentioned above), and, if it 
was anomalous for the same reasons, its K-band flux ought 
to be even more strongly suppressed, which is not seen. 



5 FUNDAMENTAL PARAMETERS FROM 
SPECTRAL MODELLING 

In this section we compare synthetic near-infrared spectra 
to the observed GNIRS spectrum of ULAS J0034-00, in 
an attempt to constrain effective temperature T e «, surface 
gravity, and metallicity. If T e g and log g can be constrained 
then mass and age can also be constrained, from evolution- 
ary models. 

The distance to ULAS J0034— 00 is currently unknown 
and therefore, in comparing against model spectra, the con- 
straint provided by the normalisation of the spectrum is 
not available. Together with the limitations of the spectral 
models (principally opacity uncertainties, described below) 
this precludes the accurate measurement of parameters T e ff, 
logg, and [m/H] for ULAS J0034— 00 by direct comparison 
of synthetic spectra to data alone. To overcome these limi- 
tations we have chosen t o use the recent detai led modelling 
of 2MASS J0415-09 by ISaumon et all |2007l ) as a bench- 
mark. Then to estimate parameters for ULAS J0034— 00, 
we compare differences in the measured spectra between 
2MASS J0415-09 and ULAS J0034-00, with differences in 
the model spectra, as parameters are changed away from the 
fiducial values for 2MASS J0415-09. Put another way, the 
models are calibrated to 2MASS J0415-09. The rationale, 
an d the method itself, are s imilar to the procedure developed 
by IBurgasser et all (|2006bl) , but with significant differences. 

We first provide details of the spectral models, and then 
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Figure 6. Direct comparison of the best-fit model with the GNIRS spectrum, where the model has been scaled to the crown of the 
J-band peak. 



present the results of the model fits. Combinations of T e g, 
logg, and [m/H] are denoted e.g. [750, 5.0, 0.0]. 



5.1 The BT-Settl models of T dwarf spectra 

The model atmospheres used to generate the theoretical 
T dwarf spectra for this comparison were generated with 
version 1 5.3 of the general-purpose stellar atmosphere code 
PHOENIX jHauschildt fc BarorJ[l999h . For the present models 
we use a setup that gives the currently best fits to observed 
spectra of M, L, and T dwarf s, updating the microphysic s 
used in the GAIA model grid (jKucinskas et al.ll2005l. l200rj) 
W e adopt the solar abundances of iGrevesse fc Sauvall 
(1998). Several recent studies of line formation in the solar 
atmosphere have determined carbon, oxyge n and nitrogen 
abundances nearly a factor of two lower (lAsplund et al. I 
12004 lAsplund et all 120051 ; lAllende Prietol 120071 ). However 
these revisions have been found to be in conflict with other 
observations and models of the solar a tmosphere, such a s 
studies of the carbon monoxide bands (|Avres et al.l [2006) . 
and in particular with helioseismology-based mod els of the 
solar interior (e.g. Dela have fc Pinsonn cault 2006|). The is- 
sue of the solar abundance scale for C, O and probably N, is 
therefore far from settled. Since the opacity in an essentially 
dust-free T dwarf atmosphere is largely dominated by the 
molecular bands of H2O, CH4, and to some extent NH3, it 
is the relative abundance of CNO that primarily determines 
the structure and spectral energy distribution of these ob- 
jects. A reduction of just these three elements by a factor of 
two is thus almost equivalen t to decreasing the overall metal- 
licity to [m/H]=— 0.3 on the lGrevesse fc Sauvall (1 19981 ) scale, 
i.e. the solar-metallicity models presented here are nearly 



identical to metal-ric h models of [m/H]=+0.3 based on the 
lAsplund et all (|2005l ) solar abundances. 

The eq uation of state (EO S) is an updated version of 
that used in lAllard et al.l l|200ll ). including ~ 10 3 species of 
atoms, ions and molecules as well as allowing for the for- 
mation and dissolution of ~ 100 grain species. One major 
update improves on the 'Cond' limit of treat ing the EOS 
and t he opacity of condensates, described in lAllard et al.l 
i|200ll ). which handled the absence of dust clouds from the 
visible photosphere at low T e ff by omitting their opacity 
from the radiative transfer calculations, while still treating 
their formation in thermodynamical and chemical equilib- 
rium. The more realistic 'Settl' treatment self-consistently 
includes the gravitational settling, sedimentation or rainout, 
of condensates, which results in a clou d layer retreating to 
larger optical depths as r e g decreases ( Allard et al. 1 120031 ; 
iHomeier et~al]|2005l ; lAllard et al. Il2007el ). In particular this 
provides a more realistic description of the depletion of re- 
fractory elements from the atmosphere. 

The water l ine op acities are taken from the calculations 
of lBarber et al.l (120061 . BT), found to be giving the best over- 
all fit to the water bands over a wide temperature range. 

Obtaining accurate and complete opacity data of CH4 
and NH3 has generally been a much greater challenge. The 
current models include the latest l inelists calculated with 
the Spherical Top Data System l|Wenger and Champion! 
1 19981 . ST DS) to account for the principal infrared bands of 
methane (|Homeier et al . 2007b). But coverage of the hot 
bands in the near-infrared, at 2.1-2.5 /im and especially at 
1.5-1.7 fj,m, is limited since the higher vibrational levels con- 
tributing to these transitions have yet to be modelled. Semi- 
empirical corrections based on high-temperature simulations 
and laboratory measurements (|Borvsow et al.l [20021 ; iBrownl 
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Figure 4. Upper: Comparison of 1-8 jim SEDs of 
ULAS J0034-00 and 2MASS J0415-09. The data for the 
latter object have been shifted by +2.5 mag. so that the SEDs 
overlap. Uncertainties are not plotted as they are comparable to 
the size of the symbols. Lower: The same information, plotted as 
the magnitude difference 7710415 — moo34 minus the magnitude 
difference in the J band, with uncertainties, showing that 
ULAS J0034— 00 is 1 mag. redder over this wavelength range. 



120051 ) have been included to account for this incomplete- 
ness. The situation is even worse for the 1.0-1.5 /im region, 
where only low-temperature linelists with mostly unknown 
ground state energies are available and opacities at brown 
dwarf temperatures can only be estimated. NH3 shows com- 
parable problems starting already at 1.9-2.1 ^m, with most 
of the lower states unknown for the 1.4-1.7 /im bands. These 
models and the quality of the opacity data are evaluated in 
more detail in comparison to the two late T dwar fs Gl 570D 
and 2MASS J0415-09 bv lHomeier et all l|2007al ). 



Figure 5. Upper: Relation between T e g and H — [4.49] colour for 
L dwarfs (filled squares) and T dwarfs (open circles) with well- 
determined temperatures. The dashed line is a linear fit to the T 
dwarfs with T cB < 1250 K. The T6p dwarf 2MASS J0937+29 is 
discussed in the text. The colour of ULAS J0034-00, H- [4.49] = 
4.00 ± 0.05, is indicated by the shaded region. Lower: Zoom of 
the upper plot, with model T e ff — colour tracks plotted for a range 
of gravity /metallicity combinations, listed in the key. The models 
are described in Section [5] The long-dashed line is the fit from 
the upper plot. 



Another important recent improvement in the PHOENIX 
models is the inclusion of new atomic line profile data based 
on more accurate interaction potentials. These profiles have 
improved the fits to observed spectra significantly over the 
standard model of a Lorentzian line shape at larg e detunings 
jAllard et al.ll2005l ; [Burrows fc Volobuvevl 12003 ). We use a 
detailed and depth-dependent line profile for each of the 
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alkali resonance lines (Li, Na, K, Rb, Cs Dl and D2, re- 
spectively) in our calculations. This is particularly impor- 
tant for the optical spectrum, being entirely dominated by 
line blanketing due to the massively broadened Nal and Kl 
doublets which in T dwarfs extend all the way out to the 
peak of the Y-band, i.e. beyond 1 /jm. In the models pre- 
se nted here we have used t he alkali line profiles describe d 



lAllard et all (12001 |2005h: lAllard fc Spiegelmanl (|2006h ; 
lAllard et al.1 i|2007bl laf) and lJohnas et al.1 (120071 ). which give 
a much improved representation of the details of these line 
shapes. 

The quality of the match of these so-called 'BT-Settl' 
models to real spectra may be gauged by examining Figs 
[2] [3] and [6] In Fig. [6] we compare the best-fit model for 
ULAS J0034— 00 determined from the spectral comparison 
below. The fit is rather good. At somewhat higher tempera- 
tures the inadequacy of the CH4 opacities manifests itself in 
the wavelength range 1.6 — 1.7 (im, as is particularly evident 
in comparing the top and bottom panels in Fig. [3] Another 
region where the fit is relatively poor is the wavelength range 
1.1 - 1.15/xm. 

In their detailed mod elling of three c ool b rown dwarfs 
ISaumon et al. I (|2006l ) and lSaumon et al.l (|2007f ) have found 
evidence for non-equilibrium effects in the carbon and ni- 
trogen chemistry of their atmospheres due to the dredge-up 
of CO and N2 from deeper and hotter layers into the pho- 
tosphere by turbulent mixing. This is parameterized by an 
eddy-diffusion coefficient K zz in their models. From analy- 
sis of the mid-infrared spectra they conclude that NH3 is 
strongly diluted in the upper atmosphere due to this effect, 
and CO is significantly enhanced, requiring K zz ~ 10 6 cm 2 
s _1 . We have included this quenching of the chemical equi- 
librium in a similar way, deriving mixing timescales from the 
same convective overshoot velocity field that is also used in 
our modelling of the cloud formation and settling processes, 
and find a similar dilution of NH3 in the mid-infrared. The 
effects of this dilution are also evident in the near-infrared 
bands of NH3, leading to considerably less absorption at 
1.45-1.55 /im and 1.9-2.1 pirn than otherwise expected. 



5.2 Spectral comparison 

To estimate the fundamental physical parameters for 
ULAS J0034— 00 we use the new model sp ectra, and follow 
a proce dure similar to that developed by iBurgasser et al.l 
(2006b]). These authors created a grid of solar-metallicity 
spectral models, for which they measured a pair of spectral 
indices: the H2O— J index (Table|3|, and an index K/H mea- 
suring the ratio of the heights of the K and H peaks. The 
same indices were measured for Gl 570D for w hich the prop- 
erties r e ff, log ,9, [M/H] are well determined (jGeballe et al.l 
120011 ; 1 aumon et al. II2006T ). The indices measured for the 
model with these properties are not identical to the indices 
measured on the actual spectrum, because the models are 
not perfect. For each of the two indices, the ratio of the ac- 
tual index over the model index provides a scaling, which 
is used to calibrate the grid of model indices. The proper- 
ties T c ff, logg, of any cool brown dwarf are then determined 
by comparing the measured indices against the calibrated 
model indices, and determining the best fit. An advantage 
of the method is that it does not require the distance of the 
object to be known, since the spectral indices record flux ra- 
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Figure 7. Plot of the normalised index Wj, quantifying the J- 
band peak width, against the normalised index K/J, measur- 
ing the ratio of the heights of the K- and J-band peaks. The 
grid shows measured indices for solar-metallicity spectral models, 
600 K < T cH < 900 K, 4.25 < logs < 5.50, and is normalised 
by associating the model [750, 5.0, 0.0] with 2MASS J0415-09, 
which is represented by the circle at (1.0, 1.0) (by definition). 
Gl 570D is represented by the square, and HD3651B by the 
hexagon. The horizontal and vertical bands plot the ±2(7 ranges 
of the normalised indices for ULAS J0034— 00. The effect of vary- 
ing metallicity is illustrated by the dashed line. 



tios. Nevertheless the effect on the measured parameters of 
varying metallicity away from solar was not explored (except 
for one source). 

Interestingly. IBurgasser et al.l l|2006bl ) estimated a value 
of T cff < 700 K for the T8 dwarf 2MASS 0939-2448 and the 
T7.5 dwarf 1114—2618, lower than the effectiv e temperature 
determ ined for 2MASS J0415-09. Recently iLeggett et all 
(2007b]) re-examined the spectra of these two sources in 
detail. They argued that the best fits were achieved with 
T c ff ~ 750 K, and sub-solar metalicity [M/H]~ —0. 3, and 
suggested that the method o f lBurgasser et al.l (|2006bT ) needs 
to be revised to take metallicity into account. 

In our compari son we also empl o y two i ndices. The first 
is the K/J index of IBurgasser et ahl ( 2006bh . defined as 



K/J = 



JF 2J 



J*. 



(1) 



where J F\ l -\ 2 denotes the integrated flux between Ai and 
A2. For the second index we chose not to use H2O — J, be- 
cause of the very limited dynamic range available at tem- 
peratures T e ff< 750 K, since the absorption is practically 
saturated. Referring to Figs[2]and[3] we define an index Wj 
which characterises the width of the J-band peak, and is 
defined by 



Wj 



J* 



(2) 



A very cool brown dwarf 13 



Table 5. Measured values of the spectral indices Wj and K/J 



Index 


Gl 570D 


2MASS J0415-09 


ULAS J0034-00 


HD 3651B 


Wj 


0.346 


0.310 


0.262±0.006 


0.342 


Wj 1 Wj(J0415) 


1.116 


1.000 


0.845±0.021 


1.103 


K/J 


0.113 


0.134 


0.128±0.003 


0.147 


K/J 1 K/J(J0415) 


0.845 


1.000 


0.960±0.022 


1.097 



where the factor of 2 in the denominator compensates for 
the wider band of the numerator. 

The measured values of the two indices for the three 
sources Gl 570D, 2MASS J0415-09, and ULAS J0034-00 
are provided in Table [S] For interest we have also pro- 
vided the measured v alues for the T7.5 dwarf HD 365 IB 
l|Mugrauer et alj|2006h . A brief discussion of this source is 
provided at the end of this section, but for the remain- 
der we consider only the first three sources. The table in- 
cludes the measured indices, as well as the values normalised 
to the measured indices of 2MASS J0415-09. The pub- 
lished spectra of Gl 570D and 2MASS J0415-09 are of 
much higher S/N than the spectrum of ULAS J0034-00 
we have obtained, so the lack of error spectra for Gl 570D 
and 2MASS J0415— 09 is unimportant in this analysis. The 
K/J indices for the three sources are all quite similar, while 
for Wj the measured index decreases through the sequence 
Gl 570D, 2MASS J0415-09, ULAS J0034-00, quantifying 
the progressive narrowing of the J-band peak visible in Figs 
and 

The normalised indices from Table [S] are plotted in 
Fig. [7] where the dotted lines outline the ±2cr ranges for 
ULAS J0034— 00, along with a grid of solar-metallicity spec- 
tral models. The grid has been scaled to the data by asso- 
ciating the model [750, 5.00, 0.0] with 2MASS J0415-09, 
and normalising the indices of all the other models to 
the measured values for this model. Choosing a different 
model to represent 2MASS J0415— 09 merely has the ef- 
fect of shifting (and slightly distorting) the grid in this 
space such that the new reference model lies at (1.0, 1.0). 
In this way it can be seen that the locations of Gl 570D and 
ULAS J0034— 00 imply temperature and gravity differences 
relative to 2MASS J0415— 09 that are, due to the regularity 
of the grid, almost independent of the actual model chosen 
to represent 2MASS J0415-09. The plot illustrates the fact 
that the two indices provide complementary information on 
the two parameters T B g and logg, at fixed metallicity. The 
effect of decreasing temperature alone is a decrease in both 
Wj and K/J, while decreasing gravity alone causes a de- 
crease in Wj, but an increase in K/J. 

If 2MASS J0415-09 and ULAS J0034-00 have simi- 
lar metallicity, the plot indicates that ULAS J0034— 00 has 
T e ff between 60 and 120 K cooler, and log g between 0.25 
and 0.5 lower. Because only two indices are used, the three 
parameters T e ff, log g, [M/H] cannot be determined inde- 
pendently. The effect of varying metallicity is illustrated in 
Fig-El by reference to the model [650, 4.75, 0.0]. The dashed 
line through this model shows the effect of varying the metal- 
licity over the range —0.2 < A[m/H]< +0.1. Because the 
line approximately follows the line of constant gravity, vary- 
ing metallicity has little effect on the temperature estimate, 



but, rather, is nearly degenerate with gravity. This degener- 
acy may be described by the relation A(log g) = — 2A[m/H]. 
Therefore ULAS J0034— 00 is of lower gravity or higher 
metallicity than 2MASS J 0415— 09 as described by the limits 
-0.5 < A(log 3 - 2[m/H]) < -0.250 

The properti e s of 2 MASS J0415-09 were determined 
by ISaumon et alj {2007) to lie along a narrow line from 
T eff , logs of 725, 5.00 to 775, 5.37. Adopting 750, 5.20 for 
2MASS J0415-09 implies 630 < T cB < 690 (±2cr range) 
for ULAS J0034— 00. The temperature estimate from the 
spectral modelling is in good agreement with the estimate 
derived from the #-[4.49] colour of 560 < T cS < 660 K 
(±2er range). Combining these two independent estimates 
gives an estimate T e ff = 648 ± 13 (la error), where the 
uncertainty reflects the random errors in the GNIRS spec- 
trum, and the mid-infrared photometry, but takes no ac- 
count of the uncertainties in the spectroscopic models, in 
the extrapolation of the temperature-colour relation (Fig. 
E>}, and in the temperature of 2MASS J0415-09. On the 
other hand since the spectroscopic models are calibrated 
to 2MASS J0415-09 at 750 K, it seems very unlikely that 
ULAS J0034-00 is hotter than 700 K. We adopt the conser- 
vative ranges 600 < T cB < 700 K, and 4.5 <logg < 5.1, as- 
suming solar metallicity, as the final estimate of the effective 
temperature and gravity of ULAS J0034— 00. Confirmation 
of this low effective temperature will require mid-infrared 
spectroscopy, and measurement of the parallax, both of 
which are scheduled. 

One of the models within the above ranges is [650, 4.5, 
0.0]. The synthetic spectrum generated with these parame- 
ters reproduces the observations quite well (given the known 
opacity problems), as shown in Fig. [S] 

Finally we note that the indices listed in Table [5] for 
Gl 570D and for HD 3651B sug gest that the latter is some 
50 K warmer. iLiu et al.l l|2007h argued that HD 3651B is 
of higher gravity than Gl 570D, in which case the relative 
values of the two indices indicate that HD 3651B must also 
have higher metallicity, in agreement with their findings. 



5.3 Distance, age and mass 

ISaumon et al.l l|2007l . their Fig. 3) plot isochrones and iso- 
mass sequences as a function of T c g and logg. They also 
demonstrate that photospheric metallicity does not signifi- 
cantly impact the evolutionary parameters in this diagram. 



6 ILiu et al.l ||2007|1 find a similar amplitude for the degeneracy be- 
tween metallicity and surface gravity determination, based on em- 
pirical comparison of the T7.5 dwarfs Gl 570D and H D 3651B and 
on exa mination of the multi-metallicity models of lBurrows et all 
ll2006h . 
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Figure 8. Comparison of the H— band spectra of ULAS J0034— 00 and 2MASS J0415— 09, at intermediate resolution. The thick solid 
line is the ISAAC spectrum of ULAS J0034— 00 binned by a factor of five, to a pixel size 4.0A. The thin solid line is the corresponding 
error spectrum. The dotted line is the spectrum of 2MASS J0415-09 from iMcLean et al.l J2003T) . 



The ranges in T eff and logg of 600-700 K and 4.5-5.1, de- 
termined above, correspond to ranges in [age (Gyr), mass 
(Afjup)] of [0.5, 15] to [8, 36]. Hence ULAS J0034-00 is likely 
to be less massive and less luminous than 2MASS J0415— 09, 
for which Saumon et al. determine a mass of 33-58 Mj up . 



The distance of ULAS J0034— 00 may be estimated from 
the modelled surface fluxes and brown dwarf radii. Normal- 
ising the models to the peak of the J band implies a dis- 
tance of 14-24 pc for the ranges in effective temperature and 
surface gravity given above. This is consistent with the ob- 
served correlation between absolute magnitud e and spectral 
type s hown for example in Figs. 8 and 9 of iKnapp et al. I 
(2004). The measured proper motion of (0'.'37 ± 0'.'07)/yr 
then translates into a tangential velocity of 25-42 kms -1 , 
suggesting m embership of the younger disk population (e.g. 
lEggen |[l99ct ). which is consistent with the age range of 0.5- 
8 Gyr derived from the temperature and gravity above. 



6 SEARCH FOR NH 3 IN THE if-BAND 

Absorption by NH3 is potentially a valuable diagnostic of 
the conditions in the a tmospheres of very cool brown dwarfs 
jBurrows et al. II2003T I. NH3 absorption was first convinc- 
ingly detected in a mid-infrared spectrum of e Indi Bab 
|Roellig et akl I2OO4T ) . and subsequent ly in all the dwarfs o f 
spectral type T2 or later, observed bv lCushing et al~l (|2006T ) 
in the mid-infrared. The search for and the interpretation 
of NH3 absorption in the near-infrared is made difficult 
by the shortcomings in the line lists referred to previously. 
ISaumon et alj (|2000T ) report the detection of very weak NH3 
absorption in the H and K bands in the spectrum of the 
T7p dwarf G1229B. The weakness of the features was inter- 
preted as a consequence of nonequilibrium chemistry. This 
finding is corroborated by analysis of NH3 absorption in th e 
mid-infrared spectrum of Gl 570D (|Saumon et al. II2OO6I ). 
At lower temperatures NH3 should strengthen, and its de- 
tection would be a valuable aid in further developing the 
spectral models of dwarfs with T h< 700 K. The NH3 ab- 
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sorpt ion may develop into a sufficiently strong spectral fea- 
ture (|Burrows et al. 1120031 ) that it will be incorporated into 
future schemes for classification beyond T8, possibly moti- 
vating tlie_CTeatioiiofa_jiew spectral type (but see discus- 
sion m iLeggett et ail (|2007bh ). Our motivation for obtaining 
a higher-resolution spectrum in the H band was to investi- 
gate whether the excess absorption seen in the low-resolution 
spectrum in the blue wing of the //-band peak is resolved 
into absorption lines of NH3 at higher resolution. 

The ISAAC spectrum of ULAS J0034-00 is reproduced 
in Fig. [8] For the purposes of the plot the spectrum has been 
binned by a factor of five, to reduce the noise, so the plotted 
spectrum is undersampled relative to critical sampling by 
a factor of about two. The thin line shows the error spec- 
trum for this binning. When further binned to match the 
sampling of the GNIRS spectrum, the two spectra are in 
close agreement, confirming the excess absorption relative 
to 2MASS J041 5-09. The dotted line in the figure plots the 
intermediate resolut ion spectrum of 2MASS J041 5— 09 from 
iMcLean et ail ||2003h . scaled to match ULAS J0034-00 over 
the wavelength range 1.56 — 1.575/xm. The overall match be- 
tween the two spectra is rather close, in terms of detailed 
features. Nevertheless, a general depression of the spectrum 
of ULAS J0034-00 shortward of 1.55/Ltm is evident. To as- 
sess the reality of any possible absorption lines we sub- 
tracted the spectrum of 2MASS J04I5-09. We then ran 
a boxcar filter of width seven (original) pixels (5.6A) over 
the difference spectrum and the variance spectrum, and cre- 
ated a S/N spectrum for an absorption line of this width. 
Three candidate absorption lines at S/N> 4 were detected, 
at wavelengths 1.5323 /im, 1.5365 /im (the most significant), 
and 1.5491 /xm. These three lines are marked in Fig. [8] Us- 
ing a broad 50 A filter, we also detect a significant feature 
centred near 1.512 /im. This lies near t he wavelength of a fea - 
ture in the NH3 laboratory spectrum (|Leggett et al.ll2~007bl ). 
although there are also many strong lines of H2O absorption 
in this region. While these lines appear to be real, and could 
be due to NH3, they are quite weak. There is no obvious 
strong feature apparent in the spectrum which would sug- 
gest that ULAS J00 34— 00 is radically di fferent, requiring a 
new spectral type. ILeggett et alj l|2007bl ) have recently ad- 
dressed the question of where the strongest NH3 features 
may appear in the spectra of very cool brown dwarfs, and 
concluded that the V-band and J-band peaks may be a more 
promising region to search. 



7 SUMMARY 

We have reported the discovery of the T dwarf 
ULAS J0034-00, discovered in UKIDSS DR1, which has 
a spectral type of T8.5 and an effective temperature of 
600 < T e fi < 700 K, making it the latest and coolest T 
dwarf known at the present time. In particular it is ~ 100 K 
cooler than 2MASS J0415— 09 and probably also has lower 
luminosity and lower mass. Despite the low inferred temper- 
ature, an intermediate-resolution spectrum in the i/-band 
failed to reveal convincingly the presence of absorption by 
NH3, although similar spectroscopy of the Y- and J-band 
peaks may be more revealing. It will be extremely valu- 
able to measure the parallax of ULAS J0034— 00, in order 
to confirm the estimated temperature without recourse to 



somewhat uncertain spectroscopic models. High-resolution 
imaging of ULAS J0034— 00 will be similarly important to 
establish whether it is isolated. 

The area covered by the UKIDSS LAS will expand by a 
factor of 20 over the next few years, so prospects for finding 
even cooler brown dwarfs are extremely good. However the 
unexpectedly blue Y — J colour of ULAS J0034-00 should 
be taken into account in designing future searches for such 
ultra-cool dwarfs. While all the current models predict that 
ultra-cool T dwarfs become bluer in J — H with decreasing 
temperature l|Hewett et al. ||2006| ). they are somewhat dis- 
crepant in their predict ions for how Y — J va ries with tem- 
perature: the models of lBurrows et al.~l (120031 ) initia lly tend 
towards redder Y — J; those of Marlev et al.l l|2002h show a 
weak trend to bluer Y — J; and the BT-Settl models show a 
stronger trend to bluer Y — J down to ~ 550 K . Given the rel- 
ative d epths in the YJHK bands in UKIDSS (|Warren et al. I 
l2007af ). searches for very cool brown dwarfs should allow for 
the possibility that they are detected only in the Y and J 
bands - one of several issues that should be clarified as more 
of these objects are discovered. 
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